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A simple model of a circularly closed dsDNA in a poor solvent is considered as an 



^ . example of a semi-flexible polymer with self-attraction. To find the ground states, 

lO \ the conformational energy is computed as a sum of the bending and torsional elastic 

components and the effective self-attraction energy. The model includes a relative 
orientation or sequence dependence of the effective attraction forces between different 
pieces of the polymer chain. Two series of conformations are analysed: a multicovered 



O 

CZ3 . circle (a toroid) and a multifold two-headed racquet. The results are presented as a 

c/3 \ diagram of state. It is suggested that the stability of particular conformations may be 

^ • controlled by proper adjustment of the primary structure. Application of the model 



to other semi-flexible polymers is considered. 
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I. INTRODUCTION 

When a double-stranded DNA (dsDNA) finds itself in a solution with a condensing agent, 
it tends to take on a conformation that would minimise its exposure to the environment and, 
correspondingly, maximise the sites of self-contact^'^. Since dsDNA is a rather stiff molecule, 
variation of its shape implies a certain energy cost which is usually described as the elastic 
energy of bending and twisting. The aim of this paper is to derive simple estimates that help 
us to decide what conformation is preferable for a circularly closed dsDNA in terms of its 
total energy. The molecule is assumed to be unknotted and torsionally relaxed when in pure 
water. We are interested in finding stable ground states and do not consider temperature 
effects. 

If we assign orientation to the centreline of the dsDNA, then two ways of condensation 
are logically possible: the remote intervals of the dsDNA chain can align next to each other 
such that their orientation is either the same or opposite. The first case is best achieved in a 
double- (or multiple-) covered circle, while the second realises as a straight line interval with 
two loops at the ends (Fig. [1]). Thus we come to familiar toroidal and racquet-like shapes^"^. 
Here we consider the series of the simplest possible conformations: a toroid and a two-headed 
racquet. The latter does not necessarily imply a significant writhe, while the first has its 
writhing number growing linearly with its number of coils^. Thus, the writhing number of 
a toroid with A^ coils equals ±(A^ ~ 1); if we neglect the thickness of the polymer, while a 
multiple-covered two-headed racquet may be made to have its writhing number vanish. We 
assume that the linking number is conserved, therefore, the non-zero writhe is compensated 
by twist which should be accounted for in the calculation of the elastic energy. Recall that, 
to avoid accumulation of twist in long ropes, sailors store them in a figure-8 fashion, a similar 
trick is used when mountaineers "butterfly" the coils of their climbing ropes. These coilings 
are analogous to the two-headed racquets. 

Change of the relative orientation of contacting pieces of dsDNA is only one, most pro- 
nounced, feature of two types of condensate. Actually, the coupling of the particular sites 
that interact with each other differ as well. On the other hand, the interaction potentials of 
two homologous and nonhomologous DNA duplexes are found to be significantly different^. 
This suggests that the monomer sequence may be intentionally tuned to amplify the effective 
attraction forces for one specific conformation. Toroids and racquets may be considered as 



analogs of the secondary structures of single-stranded nucleic acids where dsDNA plays part 
of a single stranded polymer and the hydrogen bonds between complementary nucleotides 
are replaced with the electrostatic attraction. The secondary structure of RNA or ssDNA 
may be controlled by varying their primary structure. Similarly, the primary structure may 
affect the relative stability of toroids and racquets as the dsDNA condensates. Note that the 
sequence- dependent DNA condensation was studied on double-stranded poly(dG-dC)-(dG- 
dC) (GC-DNA) and ds poly(dA-dT)-(dA-dT) (AT-DNA)I. The effects due to the relative 
orientation of the contacting pieces or to their specific sequence contents do not exist for 
such homogeneous primary structures. 

As known from earlier work (both theoretical and experimental) for open polymer 
chains^i^i^""-, the racquets can only be so-called metastable structures which means that 
they always lose to toroids though their conformational energy may not differ much from 
that of toroids. Numerical simulation suggests that interaction with adsorbing surface tend 
to stabilise metastable racquet-like structures^. Experiments with condensation of circular 
plasmids (closed un/nicked and open) also show intermediate racquet-like shapes which 
seem to occur more often for unnicked DNA—. What if the racquets be given the upper 
hand by building a DNA with specially tuned mirrored homologous subsequences in order 
to increase their interaction potential? The competing toroids must not feel it which could 
cause loss of their role as most stable structures. The model considered in this paper ex- 
plores such a possibility. To simplify analysis it will be assumed that all the contacting sites 
in racquets have uniform interaction potential which can differ from that of toroids where 
the interaction forces are also assumed to be uniform. It was recently demonstrated by 
numerical simulation that presence of a direction-dependent potential affects condensation 
of semi-fiexible polymers^. 



II. CONFORMATIONAL ENERGIES 

We use a model of the semi-fiexible polymer chain that includes two parts of the confor- 
mational energy: the elastic energy and the self-interaction energy^. The balance between 
these two components determines the geometry of the ground-state configuration. 

For a polymer molecule of length L, we describe the bending energy as that of a slender 










FIG. 1. Two ways of condesation of a circular dsDNA: A^-covered circles (toroids) (a: an initial 
conformation with N = 1, h: N = 2, c: N = 3) and M-folded two-headed racquets (d: M = 1, 
e,f,g: M = 2). Red colour marks parallel self- alignment, green antiparallel. The shapes (e), (f) 
and (g) differ with the writhing number, e: Wr = ±1, f: Wr = ±1 it 1, g: Wr = ±1 (unknotted). 



elastic rod: 



Ubend = — I K^(s)ds, 




(1) 



where B is the bending stiffness and k{s) is the curvature of the centrehne. Similarly, the 
torsional energy is computed as 
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Utors = ^ I ^\s)ds, 




(2) 



where C is the torsional stiffness and w(s) is the twist rate of the polymer. We do not 
account for dependence of the elastic moduli on the variations of the sequence of monomers. 
Neither stretching of DNA is included in the model as of minor importance. 

Each piece of the polymer chain may be in contact with up to six other pieces. Conden- 
sation is possible when every such contact decreases the conformational energy. In other 
words, the polymer prefers to contact with itself in order to screen its surface from ex- 
position to the solvent. To account for the self-attraction, we employ the concept of the 
coordination mumber un- The doubled coordination number 2a n counts the number of 
possible binding sites that are left non-occupied by A^ parallel polymer strands which fill in 
the triangular lattice in the orthogonal cross section^. There exists an exact formula for the 
minimal coordination number— 
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where [■] stands for the ceiling function: for x G M, \x\ is defined as the smallest integer 
greater than or equal to x. In particular, ai = 3, 0:2 = 5, a^ = 6, and so on. Let 7 be the 
energy gain of one contact per length. We shall assume that 7 > 0. Then for an A^-fold 
parallel bundle of length Ln the interaction energy is proportional to the difference of the 
sum of the coordination numbers of A^ separate chains and the coordination number of them 
packed together 

Un = -liaiN - aN)Lj^. (3) 

In what follows, all the lengths and energies will be presented in dimensionless form by 
normalising them on the condensation length Lc = a/-B/7 and the condensation energy 
Uc = y/B^, respectively™. We shall denote the normalised contour length of the polymer 
chain by A = L/L^. The analysis is carried out under a simplified assumption of the zero 
thickness of the polymer. 

A. Multi-covered circle (toroid) 

In this work we consider only circularly closed polymers. The reference conformation 
is the single-covered circle which minimises the elastic energy. DNA is torsionally relaxed 
so that the hnking number of its strands Ck = CUq = Twq, where Twq corresponds to the 



intrinsic twist of the B-form. The normahsed bending energy is Ubend = 27r^/A. There is no 
self- interaction and, correspondingly, Eq. (|3]) gives us zero for A^ = 1. 

The circularly closed DNA may form a A^-covered circle or a toroid without being nicked. 
The size of the toroid can only take discrete values contrary to the open case. Let p be the 
radius of the circle which is traced A^ times by the DNA before it closes onto itself so that 
L = 2iTpN. Then the curvature k = 1/ p = const and the bending energy Eq. ([1]) equals 

Ubend = -2'^U = 2tT B — . 

Since the linking number is conserved, we may estimate the change of twist as Tw — 
Twq = —Wr, where the writhing number of the A^-covered circle can be estimated as 
±(A^ — 1). We further assume that the twist is uniformly distributed along the length of the 
polymer, i.e. uj{s) = ±27t{N — 1)/L. Then the torsional energy can be found from Eq. (|2]): 
Utors = ^co'^{s)L = 27r^C ^ , and after nomalisation we have Utors = 27r^c ^ , where 
we introduced c = C/B. Note that in order to model a nicked DNA one can formally set 
c = to kill the torsional energy. 

The fact that the polymer touches itself uniformly along all its length allows us to apply 
Eq. ([3]) for Ln = L/N. Then the normalised interaction energy is simply Uint = {o^n — 
Nai)X/N. Summing up all three energies gives the total 



Ut — Ubend + ^tors + U. 
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iN^ + c{N-ir)+r^-3 X. (4) 
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The above expression coincides with Eq. (4) in Ref.— for c = 0. 

B. Two-headed racquet 

To minimise the bending, the heads should be identical. Then the elastic energy reduces 
to the doubled bending energy of a racquet head which is assumed to be a bundle of M 
identical plane inflexional elasticae^. Let the contour length of the head be denoted by x 
which is the only parameter the normalised bending energy of a single component depends 
on: Uhead = a/x, where a = A^^{ko){2kl-l), ^{ko) = 2K(A;o)-F (^^, kojA Here F denotes 
the elliptic integral of the first kind of modulus k 

z 

F(z,k)= I ^ dt, 

^ ^ J VT^^Vi - kH^ 





and K{k) = F{l,k) is the complete integral. The elliptic modulus ko is found as a root of 
the equation 2r]{ko) = ^(ko), where r]{k) = E{sn{^{k),k),k) = 2E{k) — E (^,^) and E is 
the elhptic integral of the second kind of modulus k 



E(A;) = E(l, k). Numerical solution gives ko ^ 0.8551 and a ^ 18.3331. 

To compute the interaction energy of the head, we apply Eq. (j2]) which becomes um = 
{om — Mai)x- The straight handle is a 2M-bundle of normalised length A/(2M) — x and 
Eq. ([3]) yields U2m = {c(2M — 2Ma;i)(A/(2M) — x). The total energy is the sum 2Muhead + 
Um + U2M- Like for an open polymer, the size of the racquet head may take continuous 



values. For given length, the optimal length of the head equals x = \l '^ — -a +3^/ ^^*^ 



2aM 
~OL2M- 

finally we have u = l^laMiau — (^2M + 3M) + (a;2Af/(2M) — 3)A. Note that the optimal 
racquet may only exist for polymers not shorter than IMx- Its actual size is proportional, 
for given M, to the condensation length (contrary to discrete diameters of the toroids made 
of a closed polymer). 

To account for the difference in the effective attraction forces for parallel and antipar- 
allel alignments we distinguish the interaction energy coefficient for the racquet 7,-. It is 
convenient to introduce a new parameter z/, v^ = 7^/7 which measures the orientation dif- 
ference, z/^ > 1 corresponds to stronger attraction in the antiparallel configuration (as in 
the racquet handle). On a triangluar lattice the antiparallel arrangement is frustrated and 
in the ground state one third of all neighbouring pairs have the same orientation, similarly 
to antiferromagnetic materials^. When computing the value of u, this property should be 
taken into consideration. 

To keep the common normalisation, we rewrite the racquet energy as 



Ur = 2^2aM{aM - a2A/ + 3M)z/ + (a2Af/(2M) - 3)Az/l (5) 

III. COMPARISON OF CONFORMATIONS 

Now we are able to compare the energy of the conformations as functions of the polymer 
length and the orientation parameter v"^. We begin with N- and P-covered circles which 
are members of a sequence of toroidal shapes. Comparison of the conformational energies 



(Eq. (jl])) leads to the critical length 

By use of this expression we can build a sequence of optimal toroids {Ni} with critical 
lengths {Aj}. For a polymer of length A, the A^j-circle is optimal among all toroidal shapes 
if Aj_i < A < Aj. We first set A^o = 0, Aq = and start with a circle with Ni = 1. Next, 
assuming that we know the optimal sequence up to i = / (i.e. we know Ni, i = 0, . . . ,1, and 
Aj, i = 0, . . . , / — 1) we compute A/ = XnjP* = min ^NjP, ^j+i = max P*. It is easy 

P>Nj \pf^p*=Xi 

to check that d\Ut{\,P*) < dxUt{X, Nj) for A = A/ and hence there exists e > such that 
the inequality holds for A G [A/, A/ + e). Therefore the A'^z+i-toroid is optimal for A > A/ in 
some right neighbourhood of A/. Now we can apply the same procedure to find A7+1, N1+2 
and so on. The first members of the optimal toroid sequence are (Fig. [2]): 

^t = {1,2,3,4,5, 

7,10,12,14,16, 

19,24,27,30,33, 

37,44,48,52,56,...}. (7) 

The critical lengths are shown in Fig. |3] for four different values of the torsional to bending 
stiffness ratio c (note that, as can be seen from Eq. (Q, the optimal sequence St does not 
depend on c for nonnegative values). 

The sequence Eq. ([71) may be described as a series of quintuples 

St = {3n{n + 1) + 1, 

{{3n + m + l){n + l),m= 1,2,3,4}, 
n = 0,l,2,...}. 

As proposed in Ref.— , the sequence is made with the numbers that satisfy a{N) = a{N — 
1) = a{N + 1) — 1. This characterisation must be corrected: numbers that can be expressed 
as {n + l){3n + 1), n = 1, 2, 3, . . ., must be removed from the optimal sequence. 

Next we compare the energy of the racquet-like shapes by using Eq. ([5]) for M and Q. 
We come to the equation for the critical value of the effective length of the polymer n = uX 
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FIG. 2. The sequence {Ni} for the optimal A'j-covered circles (toroids). 



for the M- and Q-racquets 
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Similarly to the toroids we build a sequence of optimal racquets {Mj} with critical lengths 
{yUj}. For a polymer of length /x, the M,-racquet is optimal among all racquet shapes if 
/ij_i < ji < Hj. To start the sequence we set Mq = 0, fiQ = 0, and Mi = 1 corresponds to 
a single-covered two-headed racquet. Then, assuming that we know the optimal sequence 
up to j = J (i.e. we know Mj, j = 0, . . . , J, and fij, j = 0, . . . , J — 1) we find /ij = 

Mj+i 



f^MjQ* = mm fiMjQ 

Q>Mj 



max Q*. As seen from Eq. ([5]), the energy of an M- 
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FIG. 3. The sequence of critical lengths {Aj} for the optimal A'j-covered circles (toroids) for four 
different values of the relative torsional stiffness c. 

racquet Ur{fi,M), for fixed u, is a linear function of fi. Moreover, it is an easy exercise 
to prove that the function M{aM — «2M + 3M) monotonically increases with M so that 
we have Ur{0,Q) > Ur{0,M) for Q > M. This means that if the graphs Ur{fi,Q) and 
Mr-(/i, M), Q > M, intersect each other at some positive fiMQ, then Ur{fi, Q) < Ur{fi, M) for 
/W > jJ'MQ- Therefore the Mj+i-racquet is optimal for /i > /ij in some right neighbourhood 
of ^j. Now we iterate the procedure to compute /ij+i, Mj+2 and so on. The first members 
of the optimal racquet sequence are: 1, 2, 3, 4, 5, 6, 7, 8, 12, 15, 18, 20, 22, 24, 26, 28, 30, 35 . . . 
(Fig. H]). The critical lengths are shown in Fig. |5l 

The last comparison to make is between an iV-covered circle and an M-racquet. Equating 
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FIG. 4. The sequence {Mj} for the optimal Mj-racquets. 



energies taken from Eqs. (jll) and (Q leads to the quadratic equation 



(3 - a2Af/(2M))/i2 - 2^/2aM{aM - a2M + 3M)/i- 
- (3 - aN/N)X'^ + 27T^{N'^ + c{N - 1)^) = 0. 



(8) 



Consider first the case A^ = 1. The coefficient of the A^ term vanishes and we have a quadratic 
equation solely for /i. It has two positive roots but only the greater one is important: 
HiM = V2{^/aM{aM - a2M + 3M) + ^/aM{aM - a2M + 3M) - 7r2(3 - a2Af/(2M)))/(3 - 
a2M/(2M)). Thus, the M-racquets win over simple circles if their length exceeds fiiM- 
Clearly, we are most interested in /in = 2(v2a + \/2a — tt^) ^ 22.4636, because for M > 2 
the simple circle is not competitive and the racquets should be compared with themselves 
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FIG. 5. The sequence of critical lengths {//j} for the optimal Mj-racquets. 

and multicovered circles. In the general case N > 2, Eq. ([8]) describes a hyperbola and it 
is convenient to resolve it for A as a function of fi. If the actual A exceeds this value, the 
energy of the M-racquet is higher than of the A^-circle. 

Our aim is to represent the results of energy comparison as a diagram of state in plane 
(/i, A). We begin with relatively short polymers (Fig. |6]). A bisector divides the first quadrant 
of the parameter plane into two triangular domains. The upper one, where A exceeds /i, 
corresponds to z/^ < 1, i.e. the stronger attraction for toroidal arrangements of DNA. Clearly, 
the racquets have no chance to compete with toroids under such a condition unfavourable 
for them. As the polymer is getting longer, the optimal condensate has a growing number of 
circle coverings. The bisector itself corresponds to the equal strength of attraction potential 
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irrespective of the DNA primary structure. When the racquet-hke shapes are benefited with 
u > 1, they appear on the surface. In any case their length should exceed the minimuni 
/ill. With growth of the polymer chain, the optimal racquets acquire more dsDNA strands. 
Similarly, toroids increase their number of coils as A becomes larger. In the {fi, A) plane 
the domain of the optimal A^j-toroid is bordered by two parallel straight lines A = A,;_i and 
A = Aj with the neighbouring Aj_i- and Aj+i-toroids. It also has a border with one or more 
racquet domains. Each piece of this border is a part of a hyperbola (Eq. ([8])). In a similar 
way, each domain of optimality of a given racquet type is bounded by two parallel lines 
/i = /ij_i and /i = fij and one or more pieces of hyperbolae (Eq. ([S])). Thus, the boundary 
between toroids and racquets is a piecewise curve made of hyperbolic intervals connecting 
triple points. The latter are either a toroid triple point where a straight line A = Aj comes 
to, or a racquet triple point where two hyperbolae meet a straight line /i = fij. The first 



triple point A (/i = /in, A = Ai = 27i\/3 + c) is special, it corresponds to the equality of 
energies of the single- and double- covered circles and the simplest racquet. 

As we have already computed the sequences {Ni}, {Mj} and the corresponding critical 
lengths, we may construct the toroid-racquet boundary by applying the following algorithm. 
Suppose we know the boundary that separates the A'/-toroids and Mj-racquets, that is a 
piece of the hyperbola that begins at a point (A/_i,/i*), /x* > /ij-i, or (A*,/ij_i), A* > Xi-i, 
depending on which coordinates are greater (both points satisfy Eq. ([8]) for N = Nj, M = 
Mj). To fix the other end of our hyperbolic piece we have to choose from two candidates: 
(1) (A/,/i**) or (2) (A*^/iJ) (again both points satisfy Eq. ([8]) for A^ = Nj, M = Mj). 
If A/ < A** (and /i** < ^j), then the next triple point is toroidal (1) and the hyperbolic 
piece ends there, otherwise the end is the racquet triple point (2). Then, to describe the 
next hyperbola piece by Eq. ([H]), we switch either to A"/+i and Mj in case (1) or to Ni and 
Mj+i in case (2). Clearly, we can iterate this procedure to build the boundary curve up to 
arbitrary length (Figs. [71 [H]). 

The values of Aj-s depend on the parameter c so that for particular values of the latter 
some triple points of different types may coincide. It looks on the diagram of state (Fig. [6]) 
like four boundary curves converge in one point B where double- and triple-toroids and single- 
and double-racquets meet all together. Actually there are a pair of triple points which are 
located indistinguishably close to each other because of the value of the parameter c chosen. 
The quadruple point only exists for a slightly different value c = 1.4830. 
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FIG. 6. The diagram of state in the plane (/x, A). The domains where the conformations shown 
are ground states are painted in different colours. The bisector u = X/n = 1 (dashed line) marks 
the case of invariant attraction forces independent of the monomer sequence. The upper triangle 
u < 1 corresponds to stronger attraction for the toroidal parallel arrangement of a polymer chain. 
The lower triangle u > 1 contains domains where racquets are ground states because their special 
geometry provides enhanced attraction. The relative torsional stiffness is fixed c = 1.5. 

As we see in Figs. [TJ [HI the boundary between the toroids and racquets remains under the 
bisector for given value of the relative torsional stiffness c. We conclude then that if there 
is no difference between the interaction energy for parallel and antiparallel alignment, then 
the racquet packing always loses to toroids. Figure [9] shows profiles of energies Ut (black) 
and Ur (white) (Eqs. (jl]) and ([5]), resp.) for i^ = 1 (the bisector at the diagram of state 
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FIG. 7. The diagram of state in the plane (/u,A) for larger length (c = 1.5). Horizontal lines 
separate the toroids and vertical the racquets. Triple points are marked on the toroid-racquet 
bondary curve. The dashed line is the bisector i' = X/fi = 1. 



in Figs. [6H8]), the regions above the energy curves are painted yellow for toroids and blue 
for racquets to clearly demonstrate the difference. Note that though the torsional stiffness 
of DNA has not been measured with satisfactory precision®, the variation of c seems to be 
within the range of ~ 0.5 to ~ 2.5. The boundary between the toroids and racquets remains 
below the bisector for such c. However, if there exists another semi-fiexible polymer with 
stronger torsional stiffness, then the racquets may have lower energy than toroids for some 
range of lengths. Indeed, as seen in Figs. [TOl [HI this occurs for c = 20, but only for a 
limited interval of lengths: A ^ 800. Still, even for such a large c, the energy of racquets is 
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FIG. 8. The diagram of state in the plane (/i, A) (c = 1.5) (a blowup of Fig. [7]). 



not significantly lower than for toroids. 



IV. DISCUSSION 

The diagram of state is in agreement with the earlier conclusion that, in the zero-thickness 
model, the racquet-like shapes can be only metastable, i.e. they can only possess slightly 
greater energy than the toroidal-like conformations. For open polymer shapes this was shown 
in Ref.-. Note that the 1-racquets (or hairpins) can turn out to be stable if the polymer (or 
a filament) has a significant thickness that cannot be ignored as for the chromatin fibre^^. 
However, as we have seen, the situation changes if the polymer interacts with itself differently 
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FIG. 9. Profiles of energies ut (black) and Ur (white) (Eqs. dH) and ([5]), resp.) for c = 1.5, u = 1. 
The regions above the energy curves are painted yellow for the toroids and blue for the racquets. 

depending on whether the sequences of the contacting base pairs are homologous or not. 
This sequence-dependent effect has been studied for DNA— . 

The circular closedness constraint affects both the toroidal and racquet conformations. 
Indeed, the racquets become less favourable because the polymer must make even number 
of end loops, one loop more than in the minimal open racquet and each loop costs extra 
bending energy. On the other hand, the non-zero writhe implies extra torsional energy for 
toroids. The latter does not happen if the DNA is nicked and its linking number is not kept 
constant. Setting zero torsional stiffness to model this case does not significantly change the 
diagram of state (Figs. [6H8]). 
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FIG. 10. Profiles of energies ut (black) and Ur (white) (Eqs. (H]) and ([5]), resp.) for c = 20.0, z^ = 1. 
The regions above the energy curves are painted yellow for the toroids and blue for the racquets. 



Another problem where a constrained semi-flexible polymer may form toroids is DNA 
under external tension in the presence of multivalent ions or polypeptides in the microma- 
nipulation experiments^'i^i"— . It was found that, depending on the enviroment, toroids can 
compete with plectonemes. 

We can estimate the condensation length for the B-form of DNA. Its bending stiffness 
B = ksTlp, where Ip ~ 50 nm. The DNA-DNA interaction energy per length can be 
estimated as 7 ~ Clfc^T nm~^-. Thus, the DNA condensation length Lc = a/-B/7 ~ 22 
nm at room temperature. For c = 1.5, the length that corresponds to the first triple point 
A is L = XiLc ~ 300 nm. Thus the shortest DNA for which we expect the stable racquet is 
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FIG. 11. Energies of the toroids ut (yellow) and the racquets Ur (blue) (Eqs. (JH) and ([5]), resp.) 
for c = 20.0, V = 1, for longer polymers. 



about 10'^ bp long. The parameter v equals ~ 1.7 which is rather high, but smaller values 
are sufficient for the racquets made of longer DNA. 

If the interacting parts happen to be homologous, then their interaction may be twice 
as strong as for nonhomologous molecules^. This effect can be used to make the racquet 
shapes stable. To this aim, the sequence may be chosen such that the homologous pieces 
come into contact in the particular racquet configuration but not in toroids. Thus we would 
have 7r = z/^7 = 27. The effective racquet condensation length then shortens by factor u 
and the normalised length /i increases by the same factor. This effect is expected to be not 
so pronounced for multicovered conformations because the sequence correlations less affect 
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the interaction energy in the tightly packed bundles^. 

As we can see on the diagram of state, there is no chance for racquets if z/^ < 1 which 
describes the situation when the attraction is enhanced for parallel orientations of the con- 
tacting dsDNA intervals. This effect could be easily achieved by building the closed DNA 
from homologous pieces joined so that their sequence orientation is same. 

On the other hand, if one of the two parts is reversed, then this would give the preference 
to the antiparallel alignment like in the handle of the 1-racquet and we have to look into 
the part of the diagram for z/^ > 1 . 

V. CONCLUDING REMARKS 

The model considered in this paper assumes that the self-interacting sites tend to arrange 
themselves in a parallel bundle. This is an idealisation that cannot be achieved in reality 
because of topological constraints--. Nevertheless in many cases DNA forms highly-organised 
hexagonal bundles^^. This justifies the parallelism approximation. 

The proposed tuning of the primary structure of dsDNA may lead to creation of intrinsi- 
cally bent pieces, e.g. the so-called A-tracts. This is known to influence the size of toroidal 
condensates^i and may also shift the balance between toroids and racquets because the com- 
peting conformations differ in curvature. The latter is constant for toroids in the infinitely 
thin approximation but varies for racquets. 

The dsDNA-dsDNA interaction potential is rather complex and it may have its minima 
for a relative orientation of the duplex axes at a non-zero skew angle which depends on 
the distance between the molecules^. This may cause formation of twisted bundles even 
for unconstrained polymers^^. In the simplest case of a single pair of interacting molecules 
the complexity of the interaction forces may lead to emergence of a plectonemic structure 
which could cause significant change of writhe. This effect can interfere with possible initial 
torsional strain of the closed DNA in its reference zero-writhe conformation^^. To release 
this strain DNA may form plectonemes which could be facilitated or hindered by complex 
self-interaction forces. Analysis of these effects is left for future work on elaboration of the 
model. 

The present model uses a rough estimate of topological features of two series of confor- 
mations focusing on the basic difference between two folding patterns. The metastability of 
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the racquet-like shapes and a possibihty to finely control their conformational energy may 
be exploited to construct a trigger mechanism which will be sensitive to small changes in 
the primary structure. 
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